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Abstract 



First-order phase transitions that proceed via nucleation of bubbles are described 
by the homogeneous nucleation theory of Langer. The nucleation rate is one of the 
{Sj \ most interesting parameters of these transitions. In previous works we have computed 

nucleation rates with two different methods: (A) employing coarse-grained potentials; 
O ', (B) by means of a saddle-point approximation, using dimensionally regularized field 

theory. In this article we compare the results of the two approaches in order to test 
their reliability and to determine the regions of applicability. We find a very good 



agreement. 



1 Introduction 



The calculation of bubble-nucleation rates during first-order phase transitions is a difficult problem with a long 
history. It has applications to a variety of physical situations, ranging from the metastability of statistical 
systems, such as supercooled liquids or vapours JjJ, to cosmological phase transitions, such as the electroweak 
or the deconfinement phase transition in QCD M. (For reviews with an extensive list of references, see refs. [0, ||, 
[| ||.) Our present understanding of the phenomenon of nucleation is based largely on the work of Langer 0. 
In the absence of impurities in the system, Langer's approach is characterized as homogeneous nucleation 
theory. The basic assumption is that the system in the metastable phase is in a state of quasi-equilibrium. 
The nucleation rate is proportional to the imaginary part of the analytic continuation of the equilibrium free 
energy into the metastable phase. The rate can be calculated by considering fluctuations of the system around 
particular configurations characterized as critical bubbles or droplets. They are saddle points of the free energy 
and drive the nucleation process. Extensive studies have been carried out within this framework in the last 
decades. (See the reviews [Q, @, [| || |5j, and rcf. j?J for work in the context of field theory.) Also alternative 
approaches have been pursued, that do not rely on the explicit introduction of droplets 

We employ here the application of Langer's approach to field theory |l(|. The basic quantity of interest is 
the nucleation rate /, which gives the probability per unit time and volume to nucleate a certain region of the 
stable phase (the true vacuum) within the metastable phase (the false vacuum) . The calculation of I relies on 
an expansion around a dominant semiclassical saddle-point that is identified with the critical bubble or droplet. 
This is a static configuration (usually assumed to be spherically symmetric) within the metastable phase whose 
interior consists of the stable phase. It has a certain radius that can be determined from the parameters of the 
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underlying theory. Bubbles slightly larger than the critical one expand rapidly, thus converting the mctastablc 
phase into the stable one. 

The bubble-nucleation rate is exponentially suppressed by the action (the free energy rescaled by the tem- 
perature) of the critical bubble. Possible deformations of the critical bubble generate a static pre-exponcntial 
factor. The leading contribution to it has the form of a ratio of fluctuation determinants and corresponds to the 
first-order correction to the semiclassical result. Apart from the static prefactor, the nucleation rate includes 
a dynamical prefactor that takes into account the expansion of bubbles after their nucleation. We concentrate 
only on the static aspects of the problem and neglect the dynamical prefactor. 

We consider a three-dimensional statistical system with one space-dependent degree of freedom described 
by a real scalar field 4>(x). For example, <p(x) may correspond to the density for the gas/liquid transition, or 
to a difference in concentrations for chemical phase transitions, or to the magnetization for the ferromagnetic 
transition. Our discussion also applies to a (3+l)-dimensional quantum field theory in quasi-thermal equilibrium 
for energy scales below the temperature. Then an effective three-dimensional description applies [Tl} |. In a 
different context our results can also be applied to the problem of quantum tunnelling in a (2+l)-dimensional 
field theory at zero temperature. 

The bubble-nucleation rate is given by J(| ^| 



Here S is the action of the system for a given configuration of the field <fi. For statistical systems the parameters 
appearing in S can be related to those in the Hamiltonian (see (l2) for an example). The action of the 
critical bubble is Sb = S (<f>b(r)) — 5(0), where 4>b{r) is the spherically-symmetric bubble configuration and 
<t> = corresponds to the false vacuum. The fluctuation determinants are evaluated either at 4> = or around 
(f> = (pb(f)- The prime in the fluctuation determinant around the bubble denotes that the three zero eigenvalues 
of the operator [5 2 S/S(j) 2 ]^ =( j >b have been removed. Their contribution generates the factor (Sb/^n) 3 ^ 2 and the 
volume factor that is absorbed in the definition of / (nucleation rate per unit volume). The quantity Eq is the 
square root of the absolute value of the unique negative eigenvalue. 

The calculation of the pre-exponential factor A is a technically difficult problem. In this letter we compare 
the results of two different methods for its evaluation in order to confirm their reliability. The first method, 
described in refs. (l^, [l^] and denoted by A in the following, is based on the notion of coarse graining and 
employs the Wilson approach to the renormalization group |l4|] in the formulation of the effective average 
action JIB)]. It can be applied to a multitude of systems described by a variety of actions. In the second method, 
described in ref. |ll| and denoted by B, the nucleation rate is calculated analytically near the so-called thin-wall 
limit. Starting from a bare action with a quartic potential, the bubble solution, its action, and the fluctuation 
determinants are obtained as power-series in an asymmetry parameter, which is small near the thin-wall limit. 
Renormalization is performed in the context of dimensional regularization near 3 dimensions. 

2 Relation between coarse graining and dimensional regularization 

In ref. JTcl nucleation rates were calculated for a model described by the bare actionF] 



This potential has the typical form relevant for first-order phase transitions in (3+l)-dimcnsional field theories 
at high temperature. Through a shift <p ~ > 4 1 + c the cubic term can be eliminated in favour of a term linear in 
(j). As a result the same potential can describe statistical systems of the Ising universality class in the presence 
of an external magnetic field. The calculation was performed in the limit that the asymmetry parameter is 
small and the two minima of the potential have nearly equal depth. The critical bubbles are not far from the 
thin-wall limit: The width of the surface is much smaller than the radius. 

1 This model has a special feature: when the two minima are taken to be exactly degenerate a Zi symmetry guarantees that they 
are equivalent. Therefore the thin- wall limit of more generic models without this particular feature can be qualitatively different. 
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We would like to verify that the predictions of method A arc in agreement with those of method B. In 
order to make the comparison we must understand how to find a coarse-grained action equivalent to a bare 
or renormalized one in the framework of dimensional regularization. This can be done by requiring that the 
two versions of the same theory describe the same physics. In the one-loop approximation we can require that 
the same effective potential is obtained in both cases. With dimensional regularization, the one-loop effective 
potential is 

Vi = V - ^-(V ") 3/2 . (3) 

Since poles in e = 3 — d only appear at the two-loop level, in the one-loop approximation the minimally 
renormalized parameters are equal to the bare ones. If desired, they can be related to physical observables via 
finite renormalization corrections jn], [HJ . 

On the other hand, the coarse graining employed in method A can be implemented by introducing an 
effective infrared cutoff that acts as a mass term ~ k 2 |l^, The coarse-grained potential 14 becomes the 
effective potential for k = 0. Only fluctuations with characteristic momenta ^ q 2 ^ k 2 contribute to the 
expression that relates Vk=o and Vk' 

v = v k=Q =v k - -^[vr /2 - (fc 2 + vnn (4) 

The one- loop effective potential V\ is complex for values of (j) such that V k n < 0. This pathology of the one-loop 
approximation could be avoided solving numerically the exact renormalization-group equation for Vk- However, 
the same problem appears at one loop in the dimensionally regularized version of the theory, eq. (^), and cancels 
out when matching the two versions of the theory Therefore, at one-loop order, the coarse-grained potential 
Vk corresponding to a given bare potential Vq employed in dimensional regularization is 

V k = V --^(k 2 + V^ 3/2 - (5) 

In the context of method A the calculation of nucleation rates is performed at values of k sufficiently large that 
the one-loop approximation for Vk, eq. (|5|), is real for all <j>. 

In method A we neglect the corrections to terms with derivatives of the field in the coarse-grained action, 
keeping only the corrections to the effective potcnatial. We have verified that including also the field-dependent 
one-loop correction to the kinetic term has a negligible influence on the following results. 



3 Comparison of the two methods 

It is convenient to use the rescalings x = x/m and <fi = 2<pm 2 /j in order to put Sq in the simplified form 



So 



4m 3 



d 3 x 



18 v 



A S - 



(6) 



We can now use h and A as parameters instead of 771,7,(7. The dimcnsionlcss parameter h = 9gm 2 /j 2 , ranging 
between and 1, controls the shape of the potential: The barrier is small for small h, while h = 1 corresponds 
to two degenerate minima. Method B is valid near the thin-wall limit, i.e. for h close to 1. The dimensionless 
parameter A = -f 2 /4m 3 = 9g/4mh controls the strength of the self-interactions of the field, and, therefore, the 
size of the loop corrections. This can be seen by rewriting eq. (||) as 



V k = V 



12tt 



(P + F ") 3/2 



(7) 



where k = k/m is a rescaled version of the coarse-graining scale k. 

Before performing the comparison we identify the range of parameters h, g where both approaches give 
reliable results. The saddle-point expansion, like any perturbative method, breaks down if the dimensionless 
coupling A is too large. In the approach A this is signalled by a strong residual dependence of the nucleation 
rate on the arbitrary coarse-graining scale k at which the calculation is performed. The origin of the problem 
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lies in the omission of the important higher-order corrections in the expansion. In the thin-wall limit, for fixed A, 
the relative importance of the pre-cxponcntial factor (that we can quantify as |lnA/ln/|, with A and / in units 
of m) is minimal and approximately amounts to a A • 20% correction. As a result, the saddle-point expansion 
is meaningful for A<2. The approximation for the saddle-point action employed in method B is valid only for 
h close to 1. As a result, the comparison of the two methods can be performed only above a certain value of 
h. This minimal value of h increases (decreases) for decreasing (increasing) values of A. The reason is that 
the prc-exponcntial factor in the nucleation rate becomes less (more) important with decreasing (increasing) A, 
so that one needs to start from a more (less) accurate approximation for the bubble action in order to have a 
meaningful comparison. For example one needs h>0.8 when A = 1. On the other hand, for very large bubbles 
(in practice when h > 0.95 i.e. when S > 10 4 ) it becomes too difficult to determine numerically the saddle point 
action and the pre-exponential factor, so that method A cannot be applied reliably. In conclusion, there is a 
range of values of h (0.8 < h < 0.95 for A = 1) and of A (A<2) where both techniques can be applied and a 
comparison is possible. 

In fig. 1 we display the bubble action (bands) and nucleation rate 
(lines) in the two approaches as a function of h for A = 1. The lower 
band depicts the value of the bubble action computed from the bare 
potential in the approach B. The band is delimited by So, its exact 
numerical value (lower thin line), and by S&, the analytical approxima- 
tion for it (upper thick line). The thickness of the band indicates the 
accuracy of the analytical approximation. The dotted line corresponds 
to the prediction of method B for the nucleation rate In/, with / in 
units of m. The upper band depicts the values of the bubble action S^k 
computed from the coarse-grained potential 14 for m 2 < fc 2 < 2m 2 . The 
fc-dependence of Sm is compensated by that of the prefactor A^k ■ As a 
result, the prediction of method A for In/, denoted by a continuous line, 
is fc-independent. The overlap of the dotted and continuous lines indi- 
cates that the predictions of the two methods for the bubble-nucleation 
rate agree when one is not too far from the thin- wall limit. 

The region of validity of both approaches leads to large bubble ac- 
tions, which makes an accurate comparison difficult in fig. 1. For this 
reason we plot in fig. 2 some useful ratios: 
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(a) (\nIy()/S^k for m 2 < k 2 < 2m 2 (thick band at the bottom); 

(b) (ln/g)/5g (continuous line at the top); 




0.5 0.6 0.7 0.8 0.9 1 
h height of the barrier) 



(c) So/Sb (dashed line); 

(d) (ln/^)/(ln/e) for m 2 < k 2 < 2m 2 (thin dark band in the middle), 

(c) (In /J prox )/(ln I B ) for m 2 < k 2 < 2m 2 (thick light band in the mid- 
dle), where j^pp iox j s obtained using the following approximation 
for the prefactor [l^] 



Figure 1 : Values of the bubble ac- 
tion and the nucleation rate in 
the two approaches. See the text 
for a description. 



mA approx 
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r 3 [K' (Mr)) -VI' (0)}dr 



1/2 



(8) 



From the form of these ratios the following conclusions can be reached: 

• Band (a) shows that the prefactor gives a significant contribution to the total nucleation rate ln/4 in 
approach A. The strong k dependence of this ratio is due to S^k ■ The rate is fc- independent to a very 
good approximation for A Sa 1. 

• Line (b) demonstrates that the prefactor is important also in approach B. 
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Figure 2: Comparison of the two methods for three different values of the coupling A. The continuous line (b) at 
the top of each plot indicates the importance of the pre] "actor in the approach B, estimated through (]hIb)/Sb- 
The thick band at the bottom (a) indicates the importance of the prefactor in the approach A, estimated as 
(In 1 \a) i 'S '_4fc and computed for m 2 < k 2 < 2m 2 . The thin dark band in the middle (d) depicts the ratio of the 
nucleation rates (m/_4)/(lnig) computed in the two approaches. The light band (e) depicts (ln/^ pplox )/(lnig), 
where In j^p tox j s obtained using the approximation in eq. (Q) . The dashed line (c) corresponds to So / Sb and 
indicates how much the approximations of method B overestimate the true saddle-point action. All dimensionful 
quantities are in units of m. 



• Line (c) indicates the region of validity of the expansion around the thin-wall approximation employed 
in approach B. The analytical expression for the bubble action Sb gives a good approximation to the 
numerically computed So only above a certain value of h that increases with decreasing A. 

• The small width of band (d) indicates that the k dependence of I11/4 is very weak. This verifies that 
approach A is reliable. The fact that band (d) is close to one, demonstrates that the two approaches agree 
very well when they are both reliable. 

• Away from the thin-wall limit, Sb overestimates the true bubble action So- At h « 0.7, for example, Sb is 
~ 10% larger than So- For A ^ 1 the approximated prefactor is small, | In Ab\ < 10% • Sb, and presumably 
has a only a ~ 10% error. As a result, method B is not sufficiently accurate and the two approaches do 
not agree. However, the fact that line (c) and band (d) deviate from 1 in exactly the same way indicates 
that the disagreement is largely caused by the inaccurate determination of So, while the estimates of the 
prefactors in the two approaches are in agreement. 

• Bands (d) and (e) are in good agreement below the value of h at which the approach B stops being reliable. 
This makes it possible to extend the region of h for which an analytical expression for the nucleation rate 
is available. Near the thin wall limit (h close to 1), the analytical expression provided by method B 
(eq. (138) or eq. (147) of rcf. |l(|) can be used. For smaller h, when the expansion around the thin- 
wall approximation breaks down (as can be checked by comparing eq. (54) of ref. [ |l6| with a numerical 
determination of the bubble action) eq. <Js|) can be used. 

Finally, we point out that in approach A fluctuations around the bubble (taken into account by the prefactor) 
enhance the nucleation rate with respect to cxp(— Sm)- On the other hand, in approach B the nucleation rate 
is smaller than its tree-level approximation, when the comparison is performed at fixed values of the bare 
parameters. This does not indicate a discrepancy. The prefactor enhances the nucleation rate in approach 
B as well, if the comparison is done at fixed values of the physical mass and couplings defined at the true 
minimum fl6|| . 



4 Conclusions 



In their joint region of validity the two methods for the calculation of the nucleation rate, based on coarse- 
grained potentials on the one hand and on dimensional regularization on the other, agree very well, thus giving 
support to their reliability. 
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